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Abstract 

We investigate the equilibrium charge distribution along a single annealed poly- 
electrolyte chain under different conditions. The coupling between the conformation 
of the chain and the local charge distribution is described for various solvent quali- 
ties and salt concentration. In salt free solution, we find a slight charge depletion in 
the central part of the chain: the charges accumulate at the ends. The effect is less 
important if salt is added to the solution since the charge inhomogeneity is localized 
close to the chain ends over a distance of order of the Debye length. In the case 
of poor solvent conditions we find a different charging of beads and strings in the 
framework of the necklace model. This inhomogeneity leads to a charge instability 
and a first order transition between spherical globules and elongated chains. 

1 Introduction 

The study of charged polymers is nowadays a field of highly active research [jjj because 
both of their relations with biological systems and of their practical applications. The 
main advantage compared to neutral polymers comes from the ability of various polyelec- 
trolytes to dissolve in polar solvents and in particular in water. The polymer backbones 
are ionized and dissociate in the solution into polyions and free counterions. It has been 
recognized that the distribution of the charges along the polyion chain is an important 
parameter which allows to distinguish two types of polyelectrolytes @, |J : the charge dis- 
tribution can be either quenched or annealed. For quenched polyelectrolytes, the charges 
along the chain are in a frozen configuration as in polystyrene sulfonate. The charge dis- 
tribution depends on the initial chemistry of the macromolecule but remains fixed in the 
whole range of pH of the solution. On the other hand, the charge distribution of polyacids 
or polybases strongly depends on pH and therefore the ionization sites are able to move 
along the chain. By tuning the pH of the solution, the total charge of the polymer can 
be monitored. One speaks of annealed or weak polyelectrolytes. This extra degree of 
freedom for the charges leads to some new and non trivial behaviors. Zhulina et al. |3J] 
have shown for grafted layers of weak polyelectrolytes the existence of a regime where the 



1 



brush thickness is decreasing when the chain grafting density is increased. This is mainly 
due to the fact that the net charge of a chain is not fixed but depends on its local envi- 
ronment. Another interesting effect is the behavior of annealed polyelectrolytes in poor 
solvent. Raphael and Joanny have predicted a first order phase transition from a compact 
conformation of the chains to an extended state ||. This phenomenon explains the sharp 
increase in the viscosity and the anomalous titration curves observed in experiments 0, |3|] . 

The purpose of the present paper is to investigate under various conditions the equi- 
librium charge distribution along one isolated annealed polyelectrolyte chain, i.e. in the 
limit of infinite dilution. This problem may seem very academic, but it is a first step 
towards a more precise understanding of collective properties such as those predicted by 
Zhulina et al.: in an annealed polyelectrolyte brush or in a related system such as a poly- 
electrolyte star, one expects the charge distribution to be non-uniform along the chain. 
This work could also be used to study related systems such as the self-assembly of weakly 
charged linear micelles, where end effects are known to have a major influence on the size 
distribution ||. 

The first estimation of the charge distribution along a polyelectrolyte chain is due to 
Berghold et al. J7| using numerical simulations. We find qualitatively the same important 
result: the charges accumulate at the chain ends where there is less neighboring charges 
than in the middle of the chain. The charge distribution is thus an end effect. However 
our model does not describe the computer simulations of ref.0 since, in their case, the 
chain conformation is almost unperturbed by the addition of electric charges. It means 
that the chains used in the simulation are not in an asymptotic regime where one can use 
the concept of electrostatic blobs. The charge distribution observed in |7j would be better 
described by a model where the neutral polymer conformation is imposed and where the 
distribution of a few charges on the chain is then investigated. We will not proceed further 
with such a model in this paper. 

This paper is organized as follows. In section we evaluate the charge distribution in a 
salt-free solution, in an ideal (theta) solvent. We consider first the effect of the finite length 
of the chain on the conformation of a polyelectrolyte with a frozen charge distribution and 
then the charge distribution along an annealed polyelectrolyte. We also generalize our 
results to the case of good solvents. In section |3|, we study the effect of added salt. Finally, 
we turn to the special case of a poor solvent, and discuss the distribution of charges in the 
so-called necklace conformation. All the results are discussed in the concluding section. 

2 Charge distribution in a salt-free solution 

2.1 Conformation of a quenched polyelectrolyte in theta solvent 

The charge distribution along a single chain is an end effect: the charges are located at 
the minima of the electrostatic potential. For an infinite chain, the average potential 
is constant along the chain; on the contrary, for a finite chain, the charges at the ends 
have less neighbors than in the center, this leads to inhomogeneities in the electrostatic 
potential. Our first step is to calculate the end effect on the conformation of a chain with a 
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frozen charge distribution. We then study the coupling between the conformation and the 
charge distribution. For a quenched chain with frozen charges, we expect a trumpetlike 
shape instead of a linear chain of constant electrostatic blobs since the chain is more 
stretched at the center than at the ends (Fig.l). 

We consider a flexible polyelectrolyte of N monomers which carries fN charges dis- 
tributed uniformly along the chain. The monomer length is denoted by b. The polymer 
backbone is described in a continuous way and the position of a monomer in space is 
represented by r(s); the first derivative |^ gives the tension along the chain. We neglect 
the effect of the counterions in this part since we are dealing with a single chain. When 
the charge distribution along the chain is annealed, the chemical potential of the charges 
related to the pH of the solution is an important variable. We present here an argument 
similar to the one given by Brochard for the shape of a tethered polymer in a strong 
solvent flow M. We write a force balance on each point of the chain which amounts to 
neglecting the fluctuations around the so-called classical path of the polymer (the most 
probable conformation). The force balance reduces to: 

where E(s) is the average electric field created by the charges of the chain. Eq. fl2.1.i|) is 
valid at every point of the chain except for the end points where the force balance on a 
section of length ds reads: 

because the chain is only stretched by one side at the end point. Since the electrostatic 
force is proportional to ds, the tension vanishes on the end point of the chain. We now 
calculate the electric field created by the polymer path r(s) in order to get a self-consistent 
equation for r(s). As we are using a continuous theory, we have to introduce short distance 
cutoffs to prevent the electric field from diverging. The relevant length at position s on 
the chain is the local stretching blob size ||: 

dz = £(z(s)) & 2 . , 

ds g(z(s)) £(z(a)) 1 " " ' 

z(s) being the polymer path projected along the direction of stretching and g(z(s)) the 
number of monomer inside the blob centered at z(s). This way of introducing the blob 
size is totally consistent with the definition of the electrostatic blob, namely the length 
below which the electrostatic interactions are not relevant. The closest charges to a test 
charge on the polymer do not contribute significantly to the electric field at this point, 
because of the natural highly fluctuating state of the neutral polymer backbone. This is 
shown more precisely in the appendix with a model of chain under tension. The linear 
charge density is |^ where Q(z) = f(^y-) 2 is the charge of the blob located at z. 

The electric field is written therefore as a sum over the contributions of all the blobs: 



E(z) = J -^- / dz 1T ^- - / d Zl —^-, (2.1.4) 
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Ib = 47r e efcT is the Bjerrum length characterizing the strength of electrostatic interactions. 
Eq.( [2.1.1|) can be integrated over the monomer position s to give an integral equation for 
the blob size: 

3kT f 2 l B kT f L/2 , ^ j f Z2 ~^ , £ 2l 
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-i/2 (.22 -Zl) Jz 2 +^ 2 

In order to obtain a simple expression for the blob size, we assume that £ 2 is slowly varying 
with z (assumption checked at posteriori), so that the blob size in the integrals is chosen 
as its value at position z. This allows for an estimate of the integrals of the right hand 
side of Eq.( FX5D : 
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where u = ls/b, and the maximal blob size £m = 
We finally obtain the blob size: 
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Where Scaling 



is electrostatic blob size one gets from the scaling theory. This 
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expression holds for values of z ranging from —L/2 + C,m to L/2 — C,m- The parameters L 
and £m are determined self-consistently by the following equations: 
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(2.1.8) 
(2.1.9) 



Eq.( |2.1.8"D is the definition of the local stretching blob size and Eq. fl2.1.9~D the definition 
of the maximal blob size. 

This functional form of the blob size is consistent with the assumption of slow varia- 
tions with z: the end effect on the conformation of the polyelectrolyte is weak. We also 
note that we recover the variational result of de Gennes et al. for the logarithmic 
correction of the total length of the chain if we neglect the z dependence of the logarithm: 
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9 ^scaling 

Although the approximation made to obtain the blob size looks quite rough, we believe 
that we catch in this way the main features of the trumpetlike shape of the chain. One 
can show that the same result could also be derived from the free energy Eq.( |A.3| ), given 
in the appendix. 
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2.2 Charge distribution on an annealed polyelectrolyte 

We now turn to the annealed charge distribution 0. The total charge of the chain is not 
fixed but it is imposed by the pH of the solution because of dissociation and recombination 
of ion pairs along the chain. The ionization sites can therefore move freely. For quenched 
polyelectrolytes, in the previous section, the charges were strictly linked to the polymer 
backbone and no motion of the charges was allowed. The only way of minimizing the 
free energy was to change the conformation. With annealed polyelectrolytes, we add an 
extra degree of freedom and the trumpet effect is less pronounced. We also expect the 
charges to accumulate preferably at the ends of the chain since the charges miss there 
approximately half of their neighbors compared to the chain center. 

In order to be more quantitative, we model the charges of the chain as a one dimen- 
sional gas of interacting particles along the polymer, in contact with a charge reservoir 
which imposes on average the local charge on a chain. The chemical potential of the 
charges is directly related to the pH of the solution |J. 

Following the results of the appendix, we write the free energy of the system: 



F [z(s) , f(s)) 
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(2.2.1) 



where f(s) is the local charge distribution along the chain, and refers to the average 
over the conformation of a chain under tension. The first two terms of Eq.( |2.2.I| ) represent 
the free energy of a quenched polyelectrolyte with an inhomogeneous charge distribution 
f(s). The third term is the entropy of an unidimensional ideal gas and the last term fixes 
the charge of the chain. The minimization of the free energy with respect to f(s) gives a 
Boltzmann law for the charge distribution: 



f(s) = /(0) exp [-(0(s) - 0(0))] 



(2.2.2) 



with 



kT 



(p(s) being the average electrostatic potential. On the other hand, the 
minimization with respect to z(s) leads to an equation for the blob size. We rewrite those 
equations in the following way: 
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(2.2.4) 
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with B = and <p{z) = (f>(z) — 0(0). The blob size £ is evaluated at the middle of the 
chain. The integrals have to be regularised by short distances cutoffs. 
Those two equations describe the coupling between the charge distribution and the con- 
formation. The reduced potential <fi(z) vanishes for an infinite chain length and must 
be small in a wide range of z values for a finite chain length; this is not true when one 
approaches the end of the chain. As a first approximation, we solve Eq. Q2.2.3| ) by iterating 
it once with <f>(z). We get a first order solution with respect to the parameter 



b 2 



[z) = In I 1 - ( — ) j (2.2.5) 



b 2 \ \L/2 

This approximation amounts to calculating the electrostatic potential of a chain of ho- 
mogeneous blobs having all the same size. The parameter /(0) is the monomer charge 
fraction at the middle of the chain, we replace it by the average charge fraction f?|: 



(f) = y[ dzf(z) (2.2.6) 

L J-L/2 

We obtain finally for the charge distribution: 




/j>) = 1 _ M/)Jscaiin g ^ ln | i _ ( _1_ ) ) +2 (l-ln2)J> (2.2.7) 
Notice that we have replaced £o by Scaling- 

Since the charge distribution f(z) depends only on the ratio 4 in this first order theory 
for the chain of blobs, one can express it easily as a function of the contour length s 
by replacing 4 by 4?. The formula Eq. fl2.2.7|) is sketched in figure 2 for different charge 
fractions. 

The expansion parameter — ^^ calm g i s the Manning parameter for the condensation of 
counterions on a chain of homogeneous electrostatic blobs (see rescaling below). In fact, 
our result is the same as the one that would obtain at fixed conformation, i.e. with a chain 
of homogeneous blobs. This is true at the first order of the theory because the end effect 
on the conformation is weak. This suggests the calculation of the charge distribution of a 
fully stretched chain, with a rodlike conformation, with the same average charge fraction. 
One finds: 

f--^H i -fe) 2 ) +2 ( 1 --)} 

We obtain therefore the important result that we can get the charge distribution of a 
chain of blobs simply by applying the following scaling or coarse-graining relations on the 
charge distribution of a fully stretched chain: 



^scaling (2.2.9) 
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(/) </).(/: (%*V (2.2.10) 



Z = *& ^ ^scaling = (^f 2 *) 1 (2-2.11) 

£ = iV& — ^scaling = Nb faf*) 1 (2-2.12) 

Those relations are valid at the first order of the theory. Up to this order, one can 
generalize this result to the case of polyelectrolyte in a good solvent in a straightforward 
way with the following relations, 



b ► Scaling (2.2.13) 

(/> — (f) 9 = (/) ^f^ 1/V (2.2.14) 
L = Nb — j U« = Nb (^f^j 1 ^ (2.2.15) 

with v — 1/2 for a gaussian chain and v — 3/5 for a chain with excluded volume. In 
this last case, £ sca iing ~ (up) 3 / 7 • ^' ie van dity of this first order theory is controlled by the 
parameter u = 1b ^^ calm B for the theta solvent case. If we use for example some typical 
parameters of weak polyelectrolytes, (/) ~ 1/20, Ib ~ b ~ 7 A , then uj = 0.4 and our 
expression should be a quite good approximation for the charge distribution. 
A better than perturbative treatment to solve Eqs. ( [2.2.31 , P-2.4]) would be required to 



obtain a more accurate expression of the charge distribution. 



3 Charge distribution with screened interactions 

We investigate in this section the effect of added monovalent salt in the solution on 
the charge distribution along the chain. The salt does not take part in the dissociation- 
recombination of the charges and the counterions: we assume therefore that the only effect 
of salt is to screen the electrostatic interactions. This may be a very rough approximation 
but it is the simplest way to take into account the salt and it allows us to get an insight 
in the charge distribution. The electrostatic potential created by an elementary charge 
reduces to: 

l B kT exp (-«r) , , 

(p(r) = 3.1 

e r 

The Debye-Hiickel screening length is given by k 2 = 8ttIbc s , c s being the concentration 
of monovalent salt. For an infinite quenched polyelectrolyte, the classical picture is a 
wormlike chain with a persistence length lp [|Kj. There are two contributions to this 
length: the bare persistence length which comes from the intrinsic stiffness of the neutral 
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backbone, and the electrostatic persistence length due to the screened interactions be- 
tween the charges that stiffen the chain. The relation between the screening length and 
this electrostatic persistence length is a controversial problem which has been discussed 
extensively for a couple of decades since the early work of Odijk [O. In this paper, we 
focus on a fully flexible chain where we have £ sca img ^ b. For this regime, the most recent 
results |T2] , [IB] seem to confirm the result of Khokhlov and Khachaturian |TJJ lp e ~ k~ 2 . 
We do not need here the precise value of the persistence length we will only assume that 

^scaling "C K 1 <C lp (3.2) 

If ^scaling ^> k -1 , the chain has a swollen conformation with an excluded volume pa- 
rameter of order k~ 2 , but in this case the inhomogeneity in the charge distribution is 
small as the electrostatic interactions are screened at length scales where they are not yet 
relevant. 

In order to derive the charge distribution, we describe the polyelectrolytes with a modified 
version of the "chain under tension model" proposed by Barrat and Joanny ||10|| . The orig- 
inal model assumes that the chain is locally stretched with a tension whose orientational 
correlations extend over distances of order lp: at short length scales, the polymer behaves 
like a chain of stretching blobs while at large length scales it shows the conformation of a 
wormlike chain. Here we focus on short distances since the charge distribution is mainly 
influenced by electrostatic interactions at these length scales. We neglect therefore the in- 
teraction between charges far apart along the chain. The main effect of those interactions 
is to induce excluded volume statistics at large length scales. 

One proceeds as before and writes the free energy for the theta solvent case: 



F[v (s)} f N ' 2 d J±fdv (s) 



kT -I -N/2 I 2& 2 \ ds 



+ f(s) (In f(s) - 1) - /W(s)l (3-3) 



Here r denotes the classical path of the polymer which follows from our simplified model 
of chain under tension. After minimization with respect to f(s), we obtain a Boltzmann 
law for the charge distribution that we rewrite it in terms of electrostatic potential only: 

±, s f^\i f N/2 a i ai i\\ /exp(-/«|r(s) - r (s')|)\ , Q A , 
<f>(s) = e,p {-) l B y^ 2 ^exp(-0( S )) { | r(3) _ r(30 | ) t (3-4) 

The use of this equation by itself without using the minimization equation with respect 
to the chain conformation allows the determination of the charge distribution at a fixed 
conformation of the polymer. We will use here the average conformation of an infinite 
chain. This approximation is justified by the considerations of section [2.1| about the end 
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effect on the conformation of a polyelectrolyte with a frozen charge distribution. Here, 
the number of charges implied in the stretching of the chain, hence the stretching force, 
are reduced because of the screened interactions. The effect on the conformation should 
be weaker than in the case of a salt-free solution. 

We evaluate the average over the chain conformation in Eq. (|3.4j) by assuming that the 
chain is stretched over distances smaller than lp by a uniform tension (kTt, t being a 
unit vector. This average is first expressed in terms of its Fourier transform: 



^ exp (-k |r(si) - r(g 2 )|) 



d 3 q 



An 



|r(sij - r{s 2 )\ U J (27r) 3 g 2 + K 2 
For distances smaller than lp, we have: 
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after averaging over the direction of stretching. 
Finally, we get |T5 
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We introduced in Eq. ( |3.7|) the following notations: ^ = 1/C is the stretching blob size. 

l2 
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z slS2 = i^{si — s 2 ) is the stretched radius associated with the contour length (s\ 



and ?% 



62 r«, 



1/2 



s 2 )j the gaussian radius of the same chain segment. The size of 
a given chain segment is the gaussian radius for distances inside a blob while at larger 
length scales, it is given by the stretched radius. For £ <C z SlS2 <C lp, we recover the 
scaling result for the chain of blobs: 



exp {—k \T{Si) 
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One can see in Eq.( |3.8| ) that the screening length is still the characteristic decay length 
of the interaction, so that two charges separated by one persistence length (lp ^> k^ 1 ) 
interact only weakly. 

According to this result and to our preceding assumption of neglecting interactions be- 
tween charges far apart along the chain, we use Eq.( |3.8| ) as an interpolating formula over 
the whole range of length scales, the short distance divergence being removed by a cutoff 
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at a size £o- in other words, since charges do not feel each other at distances larger than 
the persistence length, the charge distribution should be the same within our approxi- 
mations as the one of a chain of blobs with an imposed rodlike conformation, where the 
charges interact via a screened Coulomb potential. The conformation of the chain at 
large lengthscales is not relevant for the charge distribution. Therefore, one can compute 
easily a first order solution in (f)lB^,o/b 2 for the electrostatic potential and for the charge 
distribution: 

/W _ , + !<4^o ( Ei / / N + \ N + Ei ,,N _\_\ 2 1 



(/) b2 I VV 2 J J \\ 2 / ; Nk 

with the exponential integral Ei (x) = f +oc dt cxp (~~*) and the screening length expressed 
in terms of contour length k = ||^. The charge distribution is sketched in figure 3 for 
different screening length. We have introduce in Eq. fl3.9| ) the average charge fraction (/) 
instead of the parameter fi c ^. The characteristic decay length of the charge variation from 
the ends is of order /c -1 . 

As in the preceding section, it is interesting to compare this result to the case of a fully 
stretched chain with the same average charge : 

$- i+ W(H ta ) + *((H (3 - 10) 

This form for the charge distribution of a fully stretched chain was first derived by 
Berghold et al. 0. We obtain therefore the same physical result as in the case of a 
salt-free solution: the charge distributions Eqs.([3.9|, |3.10|) are obtained from one another 
by the following rescaling: 



b — £0 (3.11) 



(f) — (f)g = (f) [j) (3.12) 

sbn — ► - k = sbni — ] (3.13) 



9 \b 
L = Nb — ► — ^0 = Nb ^ (3.14) 



9 \b 

where ^ is the scaling electrostatic blob size. As discussed in the preceding section, those 
relations are valid both for a theta solvent (y = 1/2) and a good solvent [y = 3/5) up to 
the first order of the theory. We recover the same other qualitative conclusions as in the 
case of a salt-free solution. 

Note that our results have been derived for polyelectrolytes that fulfill the inequality 
Eq.( |3.2| ). Otherwise, one has to take more carefully into account the conformation of 
the chain. In particular, if £ ~ re -1 , one reaches the so-called cross-over regime for 
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the persistence length and a more sophisticated model is needed. We believe however 
that the characteristic length of the variation of the charge distribution will still be the 
Debye-Hiickel screening length. 

One may be tempted to interpret the simulation results of Berghold et al. in 
the framework of our theory. However, as mentioned in the introduction, they have not 
enough charges in their system to reach the asymptotic regime that we are describing. For 
example, the structure factor of the chains does not evolve when the charge fraction varies 
and is characteristic of a chain with excluded volume interactions. This indicates clearly 
that the charges do not perturb the conformation of this weakly charged polymer in a good 
solvent but that they interact anyway with each other in this highly fluctuating regime 
corresponding to a neutral polymer behavior. The charge inhomogeneity predicted by our 
model is stronger than that of the simulation. This is because their situation corresponds 
to a single electrostatic blob. The effect should be enhanced if one reaches the many blob 
regime described in this paper. 



4 Charge distribution in a poor solvent 

The distinction between quenched and annealed polyelectrolytes is particularly important 
in a poor solvent where it leads to an unusual behavior of titration curves of polyacids 
or polybases. This was first studied theoretically by Raphael and Joanny in their early 
work on annealed polyelectrolytes ||. The exchanges between the charges on the chain 
and those in the solution combined with the attractive interactions in the poor solvent 
lead to a charge instability and to a first order phase transition between spherical globules 
and elongated chains of electrostatic blobs. Throughout this section, the "electrostatic 
blob" length designs the lengthscale below which the gaussian statistics is only slightly 
perturbed by the charges. Our aim is to find the charge distribution along the chains 
when they are in the stable regime and to investigate more precisely the relationship 
between the charge distribution and the transition towards the unstable regime. In their 
work on annealed polyelectrolytes, Raphael and Joanny chose a globular blob model to 
describe the polyelectrolyte chain in a poor solvent: at short length scales the attractive 
interactions dominate over the electrostatic interaction while at larger length scales the 
electrostatic interactions dominate. Note that at very short length scales, the attractive 
interactions are not yet relevant and are dominated by the gaussian fluctuations of the 
polymer backbone. This introduces a thermal blob size £t = b/r, where r is the reduced 
temperature 1 g 1 that measures the solvent strength. The polyelectrolyte conformation 
can then be described as a chain of globular blobs, the internal structure of a globular 
blob being a close packing of thermal blobs of size £y. The monomer density in this 
collapsed state is given by r/b 3 . However, this state is always unstable with respect to 
the total length fluctuations of the chain because of the weakness of the restoring force 
to stretching of a cylindrical globule. Dobrynin et al. Wty have shown that a necklace 



conformation is more stable, where the chain can be viewed as a succession of beads 
connected by narrow strings (Fig.4). The beads have the same structure as the globular 
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blobs, while the strings are described as linear arrays of thermal blobs. We now show 
that the qualitative description of the instability of annealed polyelectrolytes with respect 
to charge fluctuations is still correct if we start from a necklace conformation instead of 
a cylindrical globule conformation. We neglect therefore in a first step all the charge 
inhomogeneities along the chain. 

The free energy of a quenched polyelectrolyte in a cylindrical globule conformation of 
length L and width D has two contributions: the surface energy = jLD, 7 being 

the surface tension of the globule given by and the electrostatic energy F^ = 
kTi B (jN) U p ^ Q a logajjthmic factor. At equilibrium, those energies are of the same order of 
magnitude and this gives the value of D ~ (up) 1 / 3 w ^ n as before u = ls/b. The cylindrical 
globule conformation does not exist at low and high values of the charge fraction /; the 
stable conformations are respectively the spherical globules and the elongated chain of 
electrostatic blobs. By comparing the free energies of those different conformations, we 
get the range of stability of the cylindrical globule: 

(^) 1/2</<tS/2m " /2 (41) 

For a quenched polyelectrolyte in a necklace conformation, we recover two analogous 
contributions to the free energy F^^ and F^\ At the level of the scaling laws, one has 

^surf ~ ^surf srnce the surface of the strings is negligible compared to that of the beads. 
The gain in free energy when one passes from the cylindrical globule to the necklace comes 
from the change in electrostatic interactions: in the necklace conformation, the beads are 
further apart. Therefore the free energy is slightly reduced for a necklace. Following the 
work of Dobrynin et a!., the equilibrium free energy reads: 



i(nec) 

= (N- M st )ru l " f/ 3 + M st -f + Nf{ru) 1 ' 2 (4.2) 
kT t 



1 /2 

In Eq.flOD, M st = fN{^)' is the total number of monomers belonging to strings. 
The first term is the self-energy of the beads. The second term is the self-energy of the 
strings. Finally, the last term represents the dominant electrostatic interactions, namely 
the interaction between beads. It is convenient to express each term of Eq. (|4.2| ) as a 
function of the ratio fe, D being the diameter of the beads equal to the globular blob size 
in the model of Dobrynin et al. This ratio is small for a typical necklace conformation 
and allows us to compare the relative magnitude of the various terms in Eq.( 



Nru l ' 3 f 2 ' 3 ~ Nt 2 ^ (4.3) 



D 



Nf{ruf/ 2 ~ ^ 2 (|) (4-4) 



5/2 

(4.5) 
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M st ^f ~ ^ 2 (|) 9/2 ( 4 - 6 ) 

Some of the terms of the free energy are not relevant, but we keep them since they help 
in the interpretation of some of the following results. 

The free energy of an annealed polyelectrolyte in a necklace conformation reads: 



p(nec) p(nec) 
r annealed _ quenched 

kT ~ kT 



+ Nf{\nf-1) -fi ch Nf (4.7) 



where we add the entropy of the charges along the chain and the chemical potential 
associated to exchanges between the charges and the solution. The equilibrium value of 
the average charge fraction / minimizes the free energy Eq.( [4.7| ): 

kT kT J + 3 N /V3 + N t J + [TU) [ ^ 

The charge chemical potential, is equal up to a constant to the pH of the solution; it fixes 
the average charge of the chain. In titration experiments, one plots fi(f) vs /. The charge 
fraction obtained from Eq. (|4.8| ) corresponds to a minimum of the free energy if 



1 2 (N — M st ) tw 1 / 3 M st u , A s 

+ 2— f->0 (4.9) 



/ 9 N / 4 /3 n r 
Up to lowest order in ^j, this inequality is rewritten as: 

/V3 > ^ rn V3 (4.10) 

Since the necklace conformation, as the cylindrical globule exists for the range of charge 
fraction given by Eq. (|4.1|) , the polyelectrolyte chain in the necklace conformation is stable 
with respect to charge fluctuations if 

r < w- 3/5 iV"- 1/5 (4.11) 

up to a numerical prefactor. One recovers the upper bound for stability of the cylindrical 
globule phase found by Raphael and Joanny. This is a trivial result since the domi- 
nant term in the free energies and so the underlying physics are the same. Notice the 
close relation between this instability and the counterion condensation on a chain with 
a quenched charge distribution. In this last case, the inequality Eq.([4.11) indicates the 



regime where one never has counterion condensation for a necklace conformation |L6] at 



fixed t. Furthermore, Dobrynin et al. |l7j show in a recent paper that there is no conden- 



sation of counterions on the beads if the inequality Eq.( |4.10| ) holds. Those analogies come 
from the fact that in our formalism for the annealed charge distribution exchanges be- 
tween the charges on the chain and in the solution are taken into account, and counterion 
condensation is also described. 
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Eq.( |4.8| ) and Eq.( |4.9| ) were derived at fixed conformation. It means that we imposed 
that M st is independent of /. If we relax this constraint and take into account the 
variation of M st with /, higher order terms are added to Eq.( |4.8|) and Eq. (|4.9Q but the 
final result is not charged. For r < u~ 3 ^ 5 A -1 / 5 , //(/) is an increasing function of /, while 
for r > w~ 3 / 5 A -1 / 5 it shows a decreasing part in a range of / values. One can show that 
the stable state for this latter case is obtained from a Maxwell equal area construction 
1 for the curve fi(f). This corresponds to two coexisting type of chains in the solution: 
spherical globules with low / and chains of thermal blobs with high /. 

A way of getting more precise information about the instability is to investigate the 
evolution of the charge distribution along the chain when the system approaches the 
transition for r < u~ 3 ^ 5 A -1 / 5 where the necklace conformation is still stable. We present 
here a simplified model for the charge distribution: we assume that each monomer of the 
beads has a charge f — Sf and each monomer of the strings has a charge / + Sf, and we 
estimate the corresponding free energy. The choice of the sign of the deviation around 
the average value / is motivated by the fact that the charges are closer in the beads than 
in the strings. If one allows the charges to move, we expect therefore a charging of the 
strings. The strength of the effect is related to the distance to the transition as shown 
below. Of course this two-phase approximation for the the charge distribution does not 
take into account the finite length of the chain as in the preceding sections but we expect 
that the effect of chain ends is small on the transition. For a sake of simplicity, we work 
as before at fixed conformation and impose M st to be independent of /. A more precise 
calculation that relaxes this constraint gives the same result at leading order. 
Since we are dealing with a homogeneous scaling theory, we use the following prescription 
for the evaluation of the changes in the free energy Eq.(fO): 

• / — ► / — Sf in the self-energy of the beads 

• /—►/ + Sf in the self-energy of the strings 

• Nf — > (A — Mat) (/ — Sf) + M st (/ + Sf) in the interaction energy between the beads 
and in the exchange term 

• Nf In / - (N - M st ) (/ - Sf) In (/ - Sf) + M st (f + Sf) In (/ + Sf) 
in the entropy of the charges 

In a first approach, if we assume that the charge fluctuations Sf and Sf are small, we 
can compute the change in the free energy AF = F^ aled (f, Sf, Sf) - F^J aled (f, 0, 0) up 
to second order in Sf and Sf and we get finally the charge distribution along the chain 
by minimizing AF with respect to Sf and Sf: 



N \ 3 /V3 z r . 

Sf = — m 4.12 

I _ 2 Til / 

/ 9 Z 4 / 3 



1 One has to equal the free energy of a spherical globule and of a chain of blobs to obtain this con- 
struction. 
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where the equilibrium charge fraction / is given by Eq. 

Before interpreting those results, we check their consistency with our assumptions. 
First the deviations y and are of order (||) with a > 1, except when the denom- 
inator of Eq. (|4.12p vanishes, a case discussed below. Therefore the effect of the charge 
distribution between the beads and the strings is weak. One can explain the weakness 
of the effect with the following argument: inside a bead, the electrostatics perturbs only 
slightly the conformation of the chain. On the opposite, two beads interact strongly 
because of their high charges. If we look for the electrostatic potential produced by a 
quenched conformation of charges, we expect the beads to behave as a potential well: the 
first relevant interactions of a test charge in a string with the neighboring charges are 
those with the two closest beads and so occur at shorter distances than for a test charge 
in a bead. In a very rough estimation, the electrostatic potential in a string is twice that 
in a bead. This explains the weakness of the charging of the strings. 
In agreement with what we were expecting, 5f and 5f are positive. The first positive 
term in the numerator of Eq. (|4.12| ) and Eq. (|4.13 ) represents the electrostatic energy of 
the charges in the beads and tends to deplete the charges from the beads while the second 
term in the numerator represents the electrostatic energy of the charges in the strings and 
tends to limit the charging of the strings. The leading term is the first one. 

As mentioned before, a change in the charge distribution of a necklace drives a change 
in the conformation. For example, if one decreases the bead charge, the bead size increases 
because more monomers carrying less charges can be aggregated. Since the charge dis- 
tribution inhomogeneities are weak, except for the region of instability, it is a reasonable 
approximation to work at fixed conformation, although we were able to manage the calcu- 
lation including a change in conformation where the parameter M st remains free to adjust 
at its optimal value as mentioned above. 

The most striking feature of Eq.( [4.12[) is that according to the inequality Eq. ([4.10 ), 
the deviation of the charge of a monomer in a bead is diverging when r reaches the 
value u~ 3//5 iV~ 1//5 . Therefore, the charge instability corresponds to a divergence of the 
inhomogeneity in the charge distribution along the chain. However, our simple calculation 
does not allow any quantitative description of the transition towards coexisting phases at 
r ~ M ~ 3 / 5 jV -1 / 5 since it is a perturbative treatment. 



5 Concluding remarks 

We have presented in this paper an estimation of the equilibrium charge distribution on an 
annealed polyelectrolyte chain under various conditions. The charge distribution is mostly 
an effect of the finite length of the chain, since an infinite chain possess a translational 
invariance (for the equilibrium properties). Our model for a salt-free solution in a theta 
or a good solvent leads to self-consistent equations coupling the conformation of the chain 
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and the charge distribution. The chain conformation is the so-called classical path of the 
polymer. All the fluctuations around this most probable path are neglected. 

When one fixes the distribution of the charges along the chain, one can investigate 
the conformation of the chain alone. We predict in this case a trumpetlike shape, with 
a local electrostatic blob size £(z) increasing towards the end of the chains. The explicit 
dependency of the blob size with z is found by using a model very similar to that used 
by Brochard to study polymers in strong flows. The end effect on the conformation is 
weak and the conformation is very close to the one obtained from homogeneous scaling 
arguments. 

When the charge distribution is annealed, we only need to consider an imposed con- 
formation for the chain as a linear array of electrostatic blobs with all the same size. The 
charges accumulate at the ends of the chain. This result does not change the qualitative 
behavior of an annealed polyelectrolyte solution compared to a solution of polyelectrolytes 
with quenched charge distributions, but it shows that annealed polyelectrolytes might 
exhibit some interesting collective properties where end effects are known to be very im- 
portant (such as adsorption on a charged surface). In the case of theta or good solvent 
for example, similar properties have been predicted [11 for weak polyelectrolyte brushes. 
With a similar model, we show that the effect of added salt is to produce a faster decay of 
the inhomogeneity in the charge distribution along the chain. This charge heterogeneity 
decays over a length of the order of the Debye-Hiickel screening length. Our predictions 
for the charge distribution cannot however be directly compared to the first simulation 
made by Berghold et al. because we are not describing the same scaling regimes for 
the charge distribution. 

In the case of a solution in a poor solvent, one finds another example of the typical 
collective properties of the annealed polyelectrolytes mentioned above: within a certain 
range of parameters, there is a first order phase transition between spherical globules 
with low / and chain of electrostatic blobs with high /, induced by a charge instability on 
the chains. We show that this instability is associated with the divergence of the charge 
distribution inhomogeneity along the chain. 

Our work suggests that it would be interesting to have more precise experimental 
or numerical evidence on the behavior of annealed polyelectrolytes in poor solvents. A 
comparison with the simulations of Micka et al. [|l£| . on the behavior of strongly charged 
quenched polyelectrolytes in poor solvent would also give an insight on the close relation 
between counterion condensation on quenched polyelectrolytes and the annealed charge 
distribution. 
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Appendix: Introduction of a blob size with a model of 
chain under tension. 



In this appendix, we show how the stretching blob size can be interpreted as a short dis- 
tance cutoff of the electrostatic interactions by considering a variational model of chain 
under tension. The description of a polyelectrolyte with such a model has proven success- 
ful (for a review, see for example and it has the important advantage of allowing for 
explicit calculations. 

The hamiltonian of a single chain is: 



kT J-n/2 I 2fe 2 \ds J 2 J-N/2 l r ( s ) ~~ r ( s ' 

We choose as a trial hamiltonian the one of a chain under a tension F(s): 



H t f N/2 , 3 fdv 



2 




kT J_ N/2 ~~ [2b 2 \ds } "•' ' ^'^ 

so that the upper bound of the free energy reads up to an irrelevant constant: 

^-C*{M# + "' — CM) 

where 




dz _ / fdr\\ F(s) F(s)b 2 



ds \ \ ds ) 1 1 F(s) 



(A.4) 



The brackets denote an average over the trial statistical weight. With Eq. flA.4 ) we choose 
z(s) as the variational parameters instead of F(s). 

To compute the electrostatic contribution to the free energy, we evaluate the average of 
the inverse of the distance between two monomers. This is done in a straightforward way: 



Erf 



/3\l/2 \z(s)-z(s')\ 
\2) \ S - S '\^ 2 b 



r(s) — r(s')\ 1 1 \z(s) — z(s')\ 



(A.5) 



with the error function Erf(x) = J* dy exp (— y 2 ). For short separations, the average 
can be rewritten as: 



r(s) -r{s')\/ t W b\s-s'\ 1/2 



(A.6) 



There is therefore no short distances divergence in the calculation of the average elec- 
trostatic potential. This is due to the highly fluctuating state of a stretched gaussian 
chain for length scales smaller than the stretching blob size. The average of the inverse 
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distances turn out to be close to the inverse of the average distances when the argument 
of the error function is of order of one. This defines in a very natural way the number of 
monomers in a blob: 




(A.7) 



Since a blob has gaussian statistics £ ~ g 1 ^ 2 b, we recover the definition of the stretching 
blob size ^ = a = 

as £ 

Therefore, we find that the charges located in a blob centered at z(s) do not introduce 
short distances divergence, and so do not contribute significantly to the average elec- 
trostatic potential at this point. For reasons of simplicity, we replace in the rest of the 
paper the regularising effect of the error function implied in the calculation of the average 
electrostatic potential by the introduction of short distances cutoffs. Finally, we notice 
that we get the same result for the blob size under the same approximations (weak effect) 
as the calculation presented in section Q by minimizing the upper bound of the free en- 
ergy derived in this appendix. We actually use this last formalism to obtain the charge 
distribution along an annealed polyelectrolyte chain. 
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Figure Captions 

Figure 1: Trumpetlike conformation: L is the total length of the chain, £m is the maximal 
blob size. 

Figure 2: Charge distribution in a salt-free solution under theta solvent conditions: N = 
130, l B = b = 7k. 

Figure 3: Charge distribution with added salt under theta solvent conditions: (/) = 1/20, 
N = 130, l B = b = 7 A. 

Figure 4: Necklace conformation: the inhomogeneity in the charge distribution is modeled 
by different charging of beads f — Sf and strings / + Sf. 
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Figure 1 Castelnovo et al. 
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Figure 3 Castelnovo et al. 



23 




Figure 4 Castelnovo et al. 



24 



